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Abstract 

We present the details of a model in General Relativity of a small 
charged black hole moving in an external gravitational and electro- 
magnetic field. The importance of our model lies in the fact that we 
can derive the equations of motion of the black hole from the Einstein- 
Maxwell vacuum field equations without encountering infinities. The 
key assumptions which we base our results upon are (a) the black 
hole is isolated and (b) near the black hole the wave fronts of the ra- 
diation generated by its motion are smoothly deformed spheres. The 
equations of motion which emerge fit the pattern of the original De 
Witt and Brehme equations of motion (after they "renormalise" ) . Our 
calculations are carried out in a coordinate system in which the null 
hypersurface histories of the wave fronts can be specified in a simple 
way, with the result that we obtain a new explicit form, particular to 
our model, for the well-known "tail term" in the equations of motion. 
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1 Introduction 



The purpose of this paper is to present a model in General Relativity of 
a small charged black hole moving in an external gravitational and elec- 
tromagnetic field. The external field is a solution of the Einstein-Maxwell 
field equations and the history of the black hole is a non-singular time-like 
world-line in this background space-time. The presence of the black hole is 
envisaged as a small perturbation of the background space-time which is sin- 
gular on this world-line (r = 0). As the world-line is approached (by letting 
r — > 0) and the mass m and charge e of the black hole are taken to be small, 
such that in the limit m — > and e — > the ratios e/r and m/r remain fi- 
nite, the perturbed gravitational field (Weyl tensor) is predominantly that of 
the Reissner-Nordstrom black hole and the perturbed electromagnetic field 
(Maxwell tensor) is predominantly the Coulomb field and the calculations 
described in this paper confirm that these limits can be achieved. These 
requirements guide us in our choice of expansions, in integer powers of r, of 
functions appearing in the metric tensor of the space-time and the potential 
1-form. 

We take the view that the assumed expansions restrict the model that 
we are constructing of a small charged black hole moving in external elec- 
tromagnetic and gravitational fields. Their generality is a topic of further 
study. Assuming that the moving small black hole is isolated (in particu- 
lar that there are no singular null geodesic generators of the histories of the 
wave fronts produced by its motion) we demonstrate how its approximate 
equations of motion are derived from the field equations as the requirement 
that the wave-fronts emerging from the moving black hole are, in the neigh- 
bourhood of the black hole, smoothly deformed 2-spheres. Approximations 
are based solely on the smallness of the mass and charge of the black hole (in 
particular slow motion is not assumed). The equations of motion derived fit 
the pattern of that of the original De Witt and Brehme [1] equations after 
they have "renormalised" . Among the significant features of our approach is 
(a) the use of a coordinate system attached to the null hypersurface histories 
of the wave-fronts of the radiation produced by the black hole motion, (b) 
the emergence of the approximate equations of motion from the requirement 
that the wave-fronts are smoothly deformed 2-spheres near the black hole, 
(c) the explicitness of the so-called tail term specific to our model and (d) 
the absence of infinities arising in our process. 

The topic of this paper has been an active area of research in general 
relativity beginning with the De Witt and Brehme [1] work (motivated by 
the classic paper by Dirac [2]). The equations of motion have been derived 
also by Beig [3] and by Barut and Villarroel For important recent work 
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on the equations of motion with radiation reaction of small black holes see 
[SI El El IE]- Some of the latter has centered on the problem of identifying 
tensor fields in the vicinity of the black hole which are singular or non- 
singular, as the case may be, on the world-line of the small black hole in 
the background space-time and the associated "regularization procedures" 
(see, for example [H]-[IH] arid the review [8]). Our work is complementary 
to the studies listed above. The precise relationship to them is a topic of 
future study however. Nevertheless recent work which should particularly 
be compared with our approach is section 4 in [5j, section 19 in [9] and [TO] , 
although the latter applies only to vacuum background geometries. Earlier 
work that has more in common with our approach, but is more specialized, 
can be found in for example [Tj5] (where the external field is considered weak) 
and in [20], [21] and [22] (where the space-time is less general than here). 

The outline of the paper is as follows: In section 2 the Einstein-Maxwell 
background space-time is described in a suitable coordinate system for our 
purposes and some consequences of the field equations required later are de- 
rived. The charged black hole space-time is introduced as a perturbation of 
the background space-time in section 3 and the consequences of imposing, 
approximately, the vacuum Einstein-Maxwell field equations are given. In 
this section the method of extracting the equations of motion of the black 
hole, using the field equations and the properties of the wave-fronts near the 
black hole, is described and the equations of motion are derived. Section 
4 is a brief discussion highlighting some properties of the equations of mo- 
tion derived in section 3. To make the paper as self-contained as possible 
calculations required for section 2 are listed in Appendix A while section 3 
requires both the calculations listed in Appendix A and the more extensive 
ones listed in Appendix B. 

2 The Background Space— Time 

We consider a small charged black hole moving in external gravitational 
and electromagnetic fields. We model the external fields by a potential 1- 
form and a space-time manifold on which it is defined which are solutions 
of the vacuum Einstein-Maxwell field equations. This space-time, which is 
otherwise unspecified in this work, contains a time-like world line (r = 0) 
on which the background Maxwell field (the Maxwell tensor field) and the 
background gravitational field (the Weyl tensor field) are non-singular. In 
the next section the small charged black hole, with mass m and charge e, 
is introduced as a perturbation of this space-time which is singular on the 
world line r = 0. The perturbed space-time will be an approximate solution 



3 



of the vacuum Einstein-Maxwell field equations having the property that in 
the limits e — > 0, m — > and r — > such that the ratios e/r and m/r remain 
finite the Maxwell field is dominated by the Coulomb field of the charge 
and the gravitational field (Weyl tensor field) is dominated by the Reissner- 
Nordstrom field of a charged black hole. To make all of these requirements 
more specific we begin by writing the line-element of the background in a 
coordinate system attached to a family of null hypersurfaces in the space- 
time having in general shear and expansion [23]. The line-element, in terms 
of a convenient basis of 1-forms, reads: 



ds* = -{&y-{$y + 2$ 6 $ A , (2.i) 



where 



d 1 = r p- 1 (e a cosh (3 dx + e~ a smhpdy + a du) , (2.2) 
d 2 = r p' 1 (e a sinh (3 dx + e~ a cosh (3 dy + bdu) , (2.3) 

tf 3 = dr + | du , (2.4) 
tf 4 = du. (2.5) 

A derivation of this line-element can be found in |24j . It is completely gen- 
eral, containing six functions p, a, j3, a, b, c of the four coordinates x, y, r, u. 
It is therefore equivalent to line-elements constructed by Sachs [25J and by 
Newman and Unti [26J. The particular form we consider here was designed 
to allow the Robinson-Trautman [27J[28J line-elements to emerge as a con- 
venient special case (this case corresponds to putting a — (3 — 0) and this 
is also a reason why the form is useful for the subject of the present pa- 
per. The hypersurfaces u = constant are null and are generated by the null 
geodesic integral curves of the vector field 8/ dr. The coordinate r is an affine 
parameter along these curves and these null geodesies have (complex) shear 

da , „ dd . « 

o = — cosh2/3 + z^ , 2.6 
or or 

and (real) expansion 

p=^-\og{rp- 1 ) . (2.7) 

We take this "background" space-time to be a solution of the vacuum 
Einstein-Maxwell field equations. We take the history of the small charged 
black hole to be a (non-singular) time-like world line in this space-time and 
we take this world-line to have equation r = 0. Assuming the line-element to 
be regular on and in the neighbourhood of r = we expand the six functions 
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introduced above in powers of r as follows (the choice of initial terms here is 
dictated by the classical work of Fermi [2H] mentioned below) : 



p = P (l + q 2 r 2 + q 3 r 3 + ...) , (2.8) 

a = a 2 r 2 + a 3 r 3 + . . . , (2.9) 

P = (3 2 r 2 + f3 3 r 3 + . . . , (2.10) 

a = di r + a 2 r 2 + . . . , (2-11) 

b = b l r + b 2 r 2 + ... ) (2.12) 

c = c + c l r + c 2 r 2 + ... . (2.13) 

In the coordinates (x, y, r, u) the potential 1-form of the background electro- 
magnetic field takes the form 

A = Ldx + Mdy + Kdu . (2.14) 

In the neighbourhood of the world line r = we expand the coefficients 
of the differentials in (12.141) in positive powers of r in such a way that the 
corresponding electromagnetic field is non-singular on r = (since this is 
the external field). A study of the exterior derivative of the 1-form f)2.14p 
reveals that the appropriate expansions of the coefficients of the differentials 
are: 

L = r 2 L 2 + r 3 L 3 + ... , (2.15) 
M = r 2 M 2 + r 3 M 3 + ... , (2.16) 
K = rK 1 + r 2 K 2 + ... . (2.17) 

The functions appearing here and in (I2.8p - fl2.13p as coefficients of the dif- 
ferent powers of r are real-valued functions of x, y, u only. We take the 
coordinates x,y,u in the range (— oo, +oo) and r in the range [0, +oo). 
Had we included in (I2.14p a term W dr we would have had to take W = 
W±r + W 2 r 2 + . . . , and this can be removed by adding a gauge term to f )2.14p 
without changing the form of the expansions fl2.15l) - fl2.17l) . 

It is easily seen from (12.61) and (12. 7p that the complex shear and expansion 
of the integral curves of d/dr are now given respectively by 

a = 2(a 2 + if3 2 )r + 3(a 3 + if3 3 )r 2 + ... , (2.18) 

p = --2q 2 r-3q 3 r 2 + ... . (2.19) 

r 

Thus near r = (for small values of r) the null hypersurfaces u = constant 
resemble future null cones with vertices on the world line r = 0. Following 
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the classical work of Fermi [2H] (see also [3U1 EI]) we can, without loss of 
generality, take the metric tensor in the neighbourhood of the world line 
r = to be the Minkowskian metric tensor, when convenient in rectangular 
Cartesian coordinates and time, up to terms of order r 2 . This has led us to 
the starting terms chosen in (12.8n — (12.13H . In addition in (12. 8p and ( 12. 13ft we 
can write (see, for example [2S] arid Appendix A) 

P = xv 1 + yv 2 + jl - \{x 2 + y 2 )) v 3 + (l + \(x 2 + y 2 )) v 4 , (2.20) 
and 
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(Q2 Q2 \ 

dx^ + dy~ 2 ) ' Cl = " 2/l0 ' ( 2 ' 21 ) 

with 

ho = = ^(logP ) , (2.22) 

-P k* = x5{ + y5i + {l-±(x 2 + y 2 )}5i-{l + ~(x 2 + y 2 )}5i. 

(2.23) 

In these formulas v l (u) is the 4- velocity of the particle with world line r = 
calculated in rectangular Cartesian coordinates and time (X 1 ) (Latin indices 
take values 1, 2, 3, 4) and u is arc length or proper-time along this line. 
rjij = diag(— 1, — 1, — 1, +1) are the components of the Minkowskian metric 
tensor in coordinates {X 1 ) and thus r)ijV l v^ = +1. Also a 1 = dv l /du is the 
4-acceleration of the particle with world-line r = (hence rjij v l a? =0) and 
on r = we have k l d/dX l = d/dr with rjijk 1 f- 7 = 1 (this latter equation 
applied to (I2.23P yields (I2.20p ). The operator A in (12.211) is the Laplacian 
on the unit 2-sphere. 

The relationship between the rectangular Cartesian coordinates and time 
(X 1 ) and the coordinates (x, y, r, u) near the world line r = is given by 
(see, for example, [32] ) 

X { = x^^ + rk* , (2.24) 

neglecting 0(r 2 )-terms. Thus the world line r = has parametric equations 
X % = x l (u) and the unit tangent vector or 4- velocity vector introduced above 
has components v % = dx l /du. We list in Appendix A some well known useful 
formulas arising from (12.241) which we will refer to in the sequel. 

The electromagnetic field calculated with the 1-form (I2.14p has the form 

F = rfA=^fA# , (2.25) 
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with F ao — —Fba and d a given by (I2.2l) - (l2.5p . Imposing Maxwell's vacuum 
field equations d*F = 0, where the star denotes the Hodge dual, we find in 
leading order (in powers of the coordinate r) the following equations to be 
satisfied by the functions L 2 , M 2 , K\\ 

2 (dL 2 dM 2 \ (dL 2 dM 2 \ 

and 

dx 2 dy { ( dx dy ) } ^ ' ^ ^ 

<9A^ 2M d( 2 fdM 2 dL 2 \) 
dy 2 dx \ \ dx dy J J 

We see immediately that the second of (12.261) is a consequence of (12.271) and 
(I2.28p . Also (12.261) implies that K\ is an I = 1 spherical harmonic: 

AiT 1 + 2K 1 = 0. (2.29) 

For our purposes the only functions of x, y, u appearing in (I2.15p -( l2.17p that 
we will require in the sequel are L 2 , M 2 and K\ and so we need not pursue 
Maxwell's vacuum field equations to higher order in powers of r. 

The leading terms in the tetrad components of the Maxwell field F a b that 
we will require are given by 
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-2 P L 2 + 0{r) , F 23 = -2 P M 2 + 0{r) , F u = K X + 0(r) . (2.30) 



On r = we can replace the basis 1-forms (I2.2I) -( J2311 by their Minkowskian 
counterparts ( 1A. 1 II) f)A. 14f) and the coordinates x, y, r, u by the coordinates 
(X*). Thus evaluating (12.301) on r = yields 



1 dk^ 1 dk j 

-F^u) k 1 —, M 2 = -F i5 {u) k l — , K 1 = Fa k % v 3 . (2.31) 



where Fij(u) = —Fjiiu) are the components of the (external) Maxwell field 
calculated in the coordinates (X*) on r = 0. We can now verify directly that 
these expressions for L 2 , M 2 , K\ satisfy the Maxwell equations (I2.26j) ( |2.29j) 
above. For example on substituting L 2 ,M 2 from (12.311) into first of (I2.26P 
yields 

Ki = p ° 2 + ^w) = l Fij{u)ki Ahj = Fl ' j{u) ¥v3 ' (2,32) 

with the last equality following from ( 1A.15I) and ( 1A.16I1 added together. 
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Taking the electromagnetic field above as source, Einstein's field equa- 
tions for the background space-time read 

Rab = 2 E ab , (2.33) 

where R ab are the components of the Ricci tensor calculated on the tetrad 
given via the 1-forms (I2.2l) -( f2~5l) . E ab are the tetrad components of the elec- 
tromagnetic energy-momentum tensor calculated using the Maxwell tensor 
(I2.25P according to the formula 

E a b = F ca F c b — - g ab F cd F cd . (2.34) 

Here g a b are the tetrad components of the metric tensor and indices on F ab 
are raised using its inverse. We will give here only the consequences of 
the field equations (I2.33P near r = which will be useful later. To satisfy 
i?33 = 2 E-33 + 0{r) we must have 

g 2 = jjp o 2 (^ + M 2 2 ) , (2.35) 

for q 2 appearing in (12. 8p . With L 2 , M 2 given by (I2.3ip . and using (IA.8I) this 
can be written 

q 2 = - 1 -FV l F V] k l ti . (2.36) 

Here F P j = F P j{u) (and F p i = r] pq F qi (u)) is the Maxwell tensor, in coordinates 
(X*) evaluated on r = 0. Now to have 

R l2 = 2E 12 + 0(r), R u - R 22 = 2 (E u - E 22 ) + O(r) , (2.37) 

requires 

2(« 2 + ^) = -|^i + *i + 4P 2 ^j , (2.38) 

where ( = x + iy and a bar will denote complex conjugation. Next, for 
A = 1 2 

R A3 = 2E A3 + 0(r) , (2.39) 

provided 

fl 1+l 6 1 + 4P 2 | = 2P^(P - 2 (a 2 + ift)) . (2.40) 
Putting (12.381) and (I2.40p together we arrive at 



>} i p^(P \a 2 + tf3 2 ))) = -(a 2 + t(3 2 ) . (2.41) 
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The approximate field equations Raa = 2 Eaa + 0(t) and P 34 = 2 E 34: + 0(r) 
both yield the same equation for the function c 2 appearing in (12.131) which will 
not be used in the sequel. The remaining field equations in the approximate 
form 

R u + Z R 24 = 2 (E 14 + iE 24 ) +0(r), Raa = 2 E M + 0(r) , (2.42) 

are automatically satisfied. To check this requires some lengthy calculations. 
With the Einstein-Maxwell vacuum field equations satisfied in the neigh- 
bourhood of the world-line r = we find the following tetrad components 
C a bcd of the Weyl conformal curvature tensor in the neighbourhood of r = 0: 



Ci3i3 + iCi323 = 6 (a 2 + ^ 2 ) + 0(r) , (2.43) 




C 3 43i+^ 3 432 = |p - 1 lai + z6i+4P^| +0(r). (2.44) 



Using the 1-forms defined in flA.lip - flA.14p and q 2 given by (I2.36P we deduce 
from these that 



6 °^' dx k dx 
1 „, „ . .• dk j , u dk 1 



a 2 = 7= P Z C W ki — kk — > ( 2 - 45 ) 



A - ^c^tf-k^, (2 ' 46) 

«! = -P 2 (<?«,■« *W — + FP.ir p . fc *^_) , (2..I.7) 




2 / dk 1 
bi = -P%[c l3kl Vvik k — + FV l F PJ k 1 ^) . (2 ..18) 

Here Cijkiiu) are the components of the Weyl tensor of this "background" 
space-time calculated on r = in the coordinates (X 1 ). Hence they satisfy 
Cijki = -C jik i = -Qjik = Ckijj and rf l C ijM = 0. Among the coefficients 
of the powers of r in fl2.8p - fl2.13p we shall only require here the functions 
92, Q!2, fa, fli, b\ given by ( 12. 36ft and fl2.45p -f j2.48l) . along with Po,co and ci 
given in (ET2D]) and (j2~2Ti) . Using the derivatives of fc* listed in (j05)l - (j07jl 



one can directly verify that the field equations ( 12.381) and (12.401) (and thus 
(J2HD) are satisfied by fl2~5oP and (J235H23EI). 



3 Black Hole Perturbation of Background 

We introduce the small charged black hole as a perturbation of the back- 
ground Einstein-Maxwell space-time described above. With the mass m 
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and charge e considered small of first order we write m = 0\ and e = 0\. 
The perturbation is introduced by modifying the expansions fl2.8p - fl2.13p as 
follows: 

p = P (l + q 2 r 2 + q 3 r 3 ... ) , (3.1) 
a = a 2 r 2 + a 3 r 3 + . . . (3.2) 
(3 = /3 2 r 2 + /3 3 r 3 + ... (3.3) 

1 a + a\ r + a 2 r 2 + . . . , (3.4) 



r 

S-i 

r 

,2 



b + b l r + b 2 r 2 + . . . , (3.5) 



e 2 2 (m + 2 /_ 



c = ^ - - v ' u ' ^ J 1 +Co + Cir + c 2 r 2 (3.6) 

The coefficients of the powers of r are functions of x, y, u. The hats are 
introduced to distinguish these coefficients from the background coefficients. 
Functions here that have non-zero background values are assumed to differ 
from their background values by Ox-terms. Functions which vanish in the 
background are assumed to have orders of magnitude: 
Oi, 6_i = 0i, bo = Oi, /-i = 02- The perturbed potential 1-form is taken 
to have the form (I2.14p with the expansions f l2.15p -( l2.17p modified to read: 

L = L + r 2 L 2 + r 3 L 3 + ... , (3.7) 
M = M + r 2 M 2 + r 3 M 3 + . . . , (3.8) 

K = (e + ^ l) + K + r K 1 + r 2 K 2 + ... . (3.9) 
r 

Again functions here which are non-zero in the background are assumed to 
differ from their background values by Oi-terms and those functions which 
vanish in the background are taken to have the orders of magnitude: L = 
Oi, M = Oi, K_i = 2 , Kq = Oi. All of these expansions have been 
chosen in order to satisfy the Reissner-Nordstrom black hole limit described 
in the Introduction. A careful study suggests that they are the minimal 
choice of expansions necessary to satisfy the black hole limit. We consider, 
however, that these assumptions restrict the model that we are constructing 
of a small black hole moving in external electromagnetic and gravitational 
fields and that their generalisation or otherwise is a topic for further study. 
In order to simplify the presentation we make the further assumption that 
near r = the potential 1-form is predominantly the Lienard-Wiechert 1- 
form (given by A = e (r _1 — ho) du up to a gauge term). This is achieved 
with the further specialisation of the functions L , M , K to satisfy 

L = O 2 , M = O 2 , K = -eh Q + O 2 , (3.10) 
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with ho given by f!2.22p . 

We emphasise that our approach is not one of "matched asymptotic ex- 
pansions" in which so-called "near zone" and "far zone" expansions of field 
variables are matched in an intermediate "buffer zone". In our perturbed 
space—time each function appearing in the metric tensor and the potential 
1-form has only one expansion in powers of r given by (3.1)-(3.6) and (3.7)- 
(3.9) and all of our deductions, including the equations of motion, are made 
using these expansions and the vacuum Einstein-Maxwell field equations. 

The expansions we have chosen here ensure that the null-hyper surf aces 
u = constant in the perturbed space-time, which will play the role of the 
histories of the wave fronts of the radiation produced by the black hole mo- 
tion, are approximately future null-cones for small r. Thus the wave fronts 
can be approximately 2-spheres near the black hole. Neglecting 0(r 4 )-terms 
the line-elements induced on these null hypersurfaces are given by 

ds 2 Q = -r 2 P - 2 (dx 2 + dy 2 ) with P = ^0(1 + ^1 + ^2 + 03) , (3.11) 

with P given by (12.201) and Qi = Oi, Q2 = O2. We shall require the per- 
turbations of these 2-spheres, described by the functions Qi and Q2 here, to 
be smooth so that as functions of x, y they are non-singular (in particular 
at x = ±00 and/or y = ±00) so that no "directional singularities" occur 
violating the notion of an isolated black hole. We will insist that in general 
in solving the Einstein-Maxwell vacuum field equations for the perturbed 
space-time and the perturbed electromagnetic field that such directional sin- 
gularities are unacceptable. Then it will follow from the Einstein-Maxwell 
field equations that necessary conditions for the 2-surfaces with line-elements 
liS.ll]) to be smooth, non-trivial deformations of 2-spheres will be the equa- 
tions of motion of the black hole. Such an origin of equations of motion can 
be found, for example, in (20j, [21] and [22] (where the space-time is less 
general than here) and in [T9] (where the external field is weak). 

The left hand sides of the vacuum Einstein-Maxwell field equations for 
the perturbed space-time are now power series in r, each with a finite num- 
ber of terms involving inverse powers of r and an infinite number of terms 
involving positive powers of r. The coefficients of the leading terms in each 
case are listed in Appendix B. They are functions of x, y, u and we require 
them to be small (in terms of e and m) with just sufficient accuracy to en- 
able us to derive the functions in the expansions fl3.ip - fl3.6l) and (13.70 - 03. 9p 
required for the equations of motion of the black hole with an 03-error (to 
include electromagnetic radiation reaction etc.). The errors that we tolerate 
to achieve this are indicated in Appendix B. Some of these errors could be 
reduced but with no influence on the equations of motion to the accuracy 
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that we calculate here. We begin by calculating the equations of motion 
with an C^-error. In this case we require from Einstein's equations (for the 
moment requiring ( IB. 331) and ( IB. 391) only to be small of second order) 

a_i = -4eP*L 2 + 2 = 1 , b-i = -4eP*M 2 + 2 = 1 , (3.12) 

and thus from Maxwell's equations (with ( IB. 211) required at this stage to be 
small of first order) we have 

^ = ^{^ + w) +0l = Kl+011 (313) 

with Ki given in (12.311) . We also need 

c = l + AQ l + 2Q 1 + 8eF ij k i v j + O 2 , (3.14) 

which emerges from the requirement that (1B.25h be small of second order. 
When these are substituted into the final Einstein field equation (given by 
requiring ( 1B.52I) for the moment to be small of second order) the differential 
equation that emerges for Qi is 

~\a(AQ 1 + 2Q 1 ) = Qma l p i -QeF ijP l v 3 +0 2 . (3.15) 

Here p % = hj k? with = <5] — v % Vj is the projection of k % orthogonal to v 1 
and satisfies p % Pi = —l,p i Vi = and since p % = k % — v % we see from (IA.15I) 
and (1A.16I) that Ap % + 2p l = so that each component of p % is an I = 1 
spherical harmonic (a spherical harmonic Q of order I is a smooth solution 
of AQ + 1(1 + 1) Q = 0). Integrating (13.151) . discarding the homogeneous 
solution with directional singularities, we have 

AQi + 2Qi =6ma i p i -6eF ij p i v j + A(u) + 2 , (3.16) 

where A(u) = 0\ is arbitrary. The first two terms on the right hand side will 
now lead to directional singularities in Q\ unless maip 1 — e p % = 2 for 
all p\ which implies the equations of motion in first approximation: 

m a>i = e Fij v 3 + 2 . (3-17) 

Hence we see the Lorentz 4-force due to the external electromagnetic field 
appearing. With (I3.17P satisfied the equation (13.161) simplifies and the di- 
rectional singularity free Q\ is a linear combination of an I = 1 spherical 
harmonic (the smooth homogeneous solution) and an / = spherical har- 
monic (the particular integral A{u)/2). In general perturbations of ( 13. lip 
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described by Q\ or Q 2 which are I = or I = 1 spherical harmonics are 
trivial. When the 2-surfaces with line-element (13.111) are viewed against the 
Euclidean 3-space background such terms in Qi or Q2 merely infinitesimally 
change the radius of the sphere (if I = 0) or infinitesimally displace the origin 
of the sphere (if I — 1) . We shall consistently discard such "perturbations" 
which means that in this case we shall take Q\ = 0. 

From now on we take Qi = in (13.111) and our objective is to solve the 
Einstein-Maxwell vacuum field equations for the perturbed space-time and 
the perturbed Maxwell field with sufficient accuracy to enable us derive the 
function Q2 in ( 13. lip and the equations of motion with an 03-error. The 
functions of x, y, u appearing in the potential 1-form (as coefficients in the 
expansions (I3.7l) - ()3.9p that are involved are L , M , K_\ and the Oi-part 
of Ki and K . Two of Maxwell's equations ( (IB. 131) and (IB.16P small of third 
order) read 



and 



These imply 




3 , (3.18) 



dy dx I V dx dy 



AK^ + 6 e 2 F tj k { v j = 3 , (3.20) 

which incidentally is ( IB. 191) . another of the Maxwell equations. This lat- 
ter equation can immediately be integrated without introducing directional 
singularities to yield 

K_i = 3 e 2 F ij k* v j + e 2 K{u) + 3 , (3.21) 

where K(u) = Oq is an arbitrary function. From (13.181) and (I3.19P we obtain 



Now the second term on the left here is an I = 1 spherical harmonic. In fact 
we can write 

1 F «l£% = *iM, (3.23) 
where F*Au) is the dual of Fij(u), on account of the relationship 

dkfe 

^— = k k —-k—, (3.24) 
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where e^ki is the Levi-Civita permutation symbol in four dimensions and k 1 
is given by f!2.23p . We integrate (13.221) to obtain 

p ° 2 (lr - If) = ~ 3 e2p ° 2 ^ k% vj + e2s{u) + ° 3 ' (3 ' 25) 

with S(u) = Oq arbitrary. This is the extent to which we shall require a 
knowledge of L and M . Substituting (I3.21j) and ( 13. 25ft into the first order 
equations ( 13.181) and ( 13.19!) verifies that they are now satisfied. We now look 
for the Oi-parts of L 2 and M%. We shall write L 2 = L 2 + h + O2 and 
M.2 = M 2 + m 2 + O2 with l 2 = 0\ and m 2 = Oi. We can obtain l 2 and m 2 
indirectly as follows: in the light of (I3.2ip and (I3.25P two of Einstein's field 
equations yield algebraically ( (1B.33h and (1B.391) required to be 3 ) 



a_i = -4eP*(L 2 + l 2 ) + 6e 2 P*F p i F pj k i — + 4e 2 P*M 2 S(u) 



: 8ti_ 

dx 

-4e 2 P 2 L 2 K{u) + 3 , (3.26) 

dk j 

6_i = -4eP 2 (M 2 + m 2 ) + 6e 2 P 2 F p i F w F — -4e 2 P 2 L 2 S(u) 

-4e 2 P 2 M 2 ir(u) +0 3 , (3.27) 
Define two functions A and P by 

a_i + 8 P^L 2 K-i = -4 e P 2 L 2 + i , (3.28) 

and 

S_i + 8 P 2 M 2 K_i = -4 e P 2 M 2 + P . (3.29) 

Now two of Einstein's field equations (( IB. 271) and ( IB. 301) both O3) can be 
written as equations for A and P neatly in the complex form 

+ iB) = -8 e 2 P 2 (L 2 + iM 2 ) 2 - 2 e 2 (a 2 + ifa) + 3 . (3.30) 

Here L 2 ,M 2 ,a 2 ,/3 2 are the background functions given by (12.311) . (12.451) and 
( I2.46P . We note that a 2 + z/? 2 can be expressed as a derivative with respect 
to ( given by ( I2.4ip . Similarly one can show that 

P 2 (L 2 + 1M2? = -2 (p o 4 (L 2 + *M 2 ) |,(L 2 + «M 2 )j . (3.31) 

Hence (13.301) can immediately be integrated giving a solution which is free 
of directional singularities: 

A + iB = 2e 2 P 4 — {P - 2 (a 2 + ^ 2 )} + 16e 2 P 4 (L 2 + ^M 2 )— (L 2 + *M 2 ) 

+e 2 U(u) + ie 2 V{u) + 3 , (3.32) 



14 



where U(u) = Oq and V(u) = Oq are functions of integration. Combining 
(PEUD with (13~2HD and (l3~27l) with (13~2^D and using the solution (gH2D means 
that we have now determined the functions l 2 and m 2 to be 

/ 2 = e F„k p ^F ij k i v1 + ±eF p i Frik i ^-±eC m k i v>k k ^ 

1 (9A- J 1 (9A-- 3 1 

+-e S(u) Fy kl ~^ + 2 e K ^ F v kl ^~A e P o 2u ( u ) + * > 

(3.33) 

m 2 = e F„^^F ij k i vi + ^eF p i F^k i ^-^eC m k i iPk k ^- 

--e Fjj fc* — + -e K{u) F tJ fc* - P ~ a V(«) + 2 . 

(3.34) 

One of Maxwell's equations (( IB. 211) required to be 2 ) now gives us directly 

K x = F i3 k* v> + P 2 + ^ - 14 e g 2 + 2 , (3.35) 

and we note that q 2 is given by (12.361) . 

We now turn our attention to the functions describing the perturbations 
of the metric tensor. We require do, bo, f-i and cq. Defining the variables 

A = d + 4 P % Kq, and B = b + 4 P 2 M 2 £T , (3.36) 

two of Einstein's equations ( (IB.28[) and (1B.31I) taken to be 2 ) can be written 
as the one complex equation 

^=(A + iB) =4mP 2 (L 2 + iM 2 ) 2 + 4m(a 2 + ip 2 ) + 2 . (3.37) 

This is integrated in similar fashion to (13.301) . The solution should have the 
property that P _1 ao and Po _1 fro are free of directional singularities since it 
is in this way that do and 60 appear in the metric tensor. The solutions of 
(I3.37P therefore that we require are 

dk j „ „„ .dkP 



a = 2eP 2 a p k p F ij k i — -2mP 2 F p i F pj k 



dx dx 
-4 m P 2 P W F v" F l3 k l — --m P^C i]kl k l v j k k — 
+X(u) + 2 , (3.38) 
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8k j 8k? 
b = 2eP 2 a p k p F lJ k i —-2mP 2 F p l F p ,k l — 

-4 m P 2 F pq k p Fij k* — - -m P 2 C m k* k k — 

+Y(u) + 2 , (3.39) 

where X(u) = 0\ and Y{u) = 0\ are functions of integration. We note that 
do and b\ possess directional singularities but P _1 a and Po^bo do not, as 
required. Another of Einstein's field equations ((B.50) required to be O3) 
provides us with the equation 

A/_x - 2e 2 a, k l + mP 2 ( ^-{p-^ x ) + ^(P^LA = 3 , (3.40) 

which using (I3.26P and (I3.27P can be rewritten as 

A(/_i - e 2 ai k i -2me F^ k { v j ) = 3 , (3.41) 

from which we conclude that 

/_! = e 2 ai k i + 2me Fy k { v j + G(u) + 3 , (3.42) 

with G(u) = 2 a function of integration. Now we obtain directly from an 
Einstein equation ( (1B.25|) taken 3 ) the function Cq. After some considerable 
simplification it can be written 

c = 1 + 8 e k { v 3 + AQ 2 + 2 Q 2 - 4 e 2 C m k l v 3 k k v l 
+8 e 2 F ij F^ + 24 e 2 ^- k* v j ) 2 - 24 e 2 F p ; F pj k l v j 

65 d 
+— e 2 F P i F pj k l k ] + 16 e 2 K(u) F %j k l v J + 4 e 2 U{u) ~-(log P ) 

+4eV( M )|-(logPo) + 03 . (3.43) 



dy 

The perturbation in the function c± = c\ + 0\ is now obtained from ( IB. 261) 
or ( IB. 491) required to be 2 . We shall not need this function in the sequel. 
The remaining quantities in (IB. 13[) (1B.21I) and (1B.23j) - (1B.51j) can now be 
evaluated and their orders of magnitude are indicated in Appendix B. The 
differential equation for Q 2 emerges from the remaining quantity (lB~52l) re- 
quired to be O3. This expression is lengthy and so we split it into a sum of 
three terms which we call T\,T 2 ,T 3 . The first of these terms involves only 
d_i and b-\ and is given by 

_ ( d *-A 2 ( dL A 2 _ I ( d *-A 2 _ I (m-A 2 

1 I dx J I dy J 2 \ dy I 2 I dx ) 
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<9y cte <9x By Bx <% 

+ 2L > w P ^ ~ 2 (a - l)2p ° ^ <p ' l) - 2 (Ll)2p ° 

dx dy \ dx 2 dxdy J 

To evaluate this we only require the leading terms in d__i and 6„i given in 
fl3~26l) and fl3~27l) . By fl230l these can be written 

a_i = 2 e Pq^is + 2 , and Li = 2 e P F 23 + 2 . (3.45) 

Direct calculation reveals 

{BP BP 1 

F 13 - F 23 + F 12 1 + 2 , (3.46) 

2 = 2e |^ F 13 + ^ F 23 - F 34 j + 2 , (3.47) 

(3.49) 

where _F a b are the components of the external electromagnetic field evaluated 
on r = and referred to the tetrad given by flA.lip - flA.14p . Substitution in 
f)3.44p leads to the simplification 

T x = 4 e 2 {(F 13 ) 2 + (F 23 ) 2 } + 8 e 2 {F 13 F 14 + F 23 F 24 } - 8 e 2 (F 34 ) 2 + 3 . 

(3.50) 



(9a_i 


db^x 


By 


Bx 


Ba-i 


db-! ^ 


Bx 


By 


9 2 a_i 




<9x By 


■{- 



p 

-^34 , 
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This is further reduced by noting that 

(Pis) 2 + (F 23 ) 2 = -F\ F pj k* k j , (3.51) 

Pis F 14 + P 23 P 24 = —F p i F pj k l + ~ F\ F pj V W - {F tJ k l v^f , (3.52) 
and 

F M = F ij k i v i . (3.53) 

Hence we can write 

T x = -8 e 2 F P i F pj k l v j - 16 e 2 (F ij k i v j ) 2 . (3.54) 

The second term T 2 is chosen because it involves the function cq explicitly. 
It is defined by 

(3.55) 

Explicit evaluation yields 

T 2 = -6^^^ + ^A(Ag 2 + 2g 2 )-3e 2 t/(n)^(logPo) 

-3 e 2 V(u) J?- (log P ) - 12 e 2 P:(w) F„ k* v j + 11 e 2 C ijfc; fc* w J A; fc 



59 



- — e 2 P ij - 42 e 2 (P ij fc* ^') 2 - 48 e 2 F p { F PJ 



v v 



-59 e 2 F p i F pj k l y + — e 2 F P i F pj k l v j + 3 . (3.56) 



The final term T 3 consists of the terms remaining in (1B.52[) . It is given by 
T 3 = 3mP 2 |^(P - 2 a ) + ^(P - 2 6 )| -4^ + 12/^ /lo + em/io 

,2 r 5i ^o , Q _ d2 dK A D2r , d2 ,, dK_ x 



+8mP z L 2 —^ + 8m P Z M 2 — - - 4 P Z L 2 —^- - 4 P 2 M 2 



^ A -dK X dK_ X 



dx dx <9y <9y 



\e 2 PZli-(Po 2 ai) + ^.(Po 2 bi)) • (3-57) 
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When this is evaluated with the functions derived above it results in (using 
h = di k l again for convenience) 



6 m h + 2 e 2 ai a { + 14 e 2 h 2 - 4 e 2 h - 4 G + 12 h G - 8 m e F tj k l v j 
+ (e m 2 + - e 2 ) F lj Fij - (12 m 2 + 5 e 2 ) C ijfci F k k v l 

d d 
-6mX(u) — (logP ) - 6mY(u) — (logP ) 

+ (50 e 2 - 36 m 2 ) (Fy U v j f + 12 e 2 F P i F t 
14 



(3.58) 



p.i 



V V 



36 m' L 

+o 3 . 



— e 2 F p i F pj k l v 1 + (18 m 2 + 5 e 2 ) F P i F pj k l W 



The approximate field equation giving us a differential equation for Q 2 is 
obtained from 7\ + T 2 + T 3 = 3 with Ti,T 2 , T 3 given by f)3.54p . (13.561) and 
(I3.58p . Making use of the unit space-like vector field p l introduced prior to 
(I3.16P we can write the result in the form 



where 



-A(AQ 2 + 2Q 2 ) = A Q + A, + A 2 + 3 , 



16 

A = -4G--e 2 F p i F pj v i v j 



(3.59) 
(3.60) 



which is an / = spherical harmonic, 



A\ = Qmaip 1 — QeF^p 1 v j — Ae 2 h k a k p l — 8 me F^p 1 v j 
-8 e 2 h\ F p k F pj p i v j + 12 G a { p { - 12 e 2 K{u) p { v j 

-3 (e 2 U(u) + 2mX(u)) ^ogP ) 
-3 (e 2 V(u) + 2mY(u)) -^(logPo) , 



(3.61) 



which is an I = 1 spherical harmonic, and 

A 2 = 6 (e 2 - 2m 2 ) - 4 (2 e 2 + 9m 2 ) X(2) + 18(m 2 - 3e 3 ) X(3) + 18e 2 x (4 ) , 

(3.62) 

which is an / = 2 spherical harmonic. Any / = 2 spherical harmonic can be 
written 

X = Dij(u) k { k j , (3.63) 

with 



Da = D 







rf j D i 







(3.64) 



19 



In (I3.62p each X{a) for a = 1, 2, 3, 4 has the form 

X( a ) = (a)Dij k l W , (3.65) 

with 

(1) Aj = Ojyi^v 1 , (3.66) 

(2) A, - W-^^^', (3-67) 

(3) Aj = F^F^h^h)- 1 -^^ 1 F%F pl , (3.68) 

(4) Aj = a,i cij - ^h i:j a k a h . (3.69) 

The last two terms in (I3.6ip can be simplified by noting that from (I2.20p and 
( I2.23P we can write 

^(logP ) =c i k i = c iV l , ^( lQ g p o) = = d lP l , (3.70) 

with 

(k= (^(v 3 -v^,0,~v\ V) and efc = (o, ^(i; 3 - v% -^ 2 , ^ 2 

(3.71) 

covariant vectors (in coordinates {X 1 )) defined along r = and everywhere 
orthogonal to v l . Provided A = 3 , (I3.59P can be integrated without the 
introduction of directional singularities to read 

AQ 2 +2Q 2 = A 1 + ^A 2 + 3 , (3.72) 

up to the addition of an arbitrary function of u which makes a trivial con- 
tribution to Q2- For Q 2 to be free of directional singularities we must have 
A\ = 3 for all p % such that p l Vi = and this leads to the equations of 
motion 



mat = e F {j v j + -e 2 h\d k + -e 2 h\ F p k F pj v j 

4 „ 
+ -m e Fij v J -2Gdi + 2 e 2 K(u) F y v 3 + U (u) a + V(u) di + O s , 

(3.73) 

where we have written U = (e 2 U(u) + 2mX(u))/2 = O2 and V = (e 2 V(u) + 
2mY(u))/2 = O2. Putting A\ = O3 and discarding the geometrically trivial 
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homogeneous solution of (13.721) we see that Q 2 = —A 2 /12 and this describes 
smooth non-trivial perturbations of the wave fronts near the black hole as 
required. G{u) is given by (13.601) with A = O3 and thus can be writ- 
ten as a definite integral for which we would naturally choose (— 00, u] as 
the range of integration resulting in —2 G cii in (13.731) contributing to a 'tail 
term'. We can write U Ci + V d{ = QijV^ with f^- = — vanishing except 
for Q 13 = Q 41 — U/2 and Q 23 = ^42 = V/2. Define uiij(u) = —Uji{u) by 
Uij = 2e 2 K(u)Fij + Qij. If we now make the 1 -parameter family of in- 
finitesimal Lorentz transformations on the unit tangent vector, v 1 — > v l = 
v l — (4/3)m e F*.,- f J — (l/m)uj l j f J , we find that, after dropping the bar on v 
and its derivatives with respect to u, the equations of motion take the form 

m a t = e F tj v j + ^e 2 h\a k + ^e 2 ^ F p k F pj v i + % + 0^ (3.74) 

where 

% = - {u k j Fji - Ui j F jk - 2 G F ik ] v k = 2 . (3.75) 
m 

Since cjy can be written as an integral, whose range we would naturally take 
to be (—00, it] as in the case of G above, we see that (13. 75ft is a 'tail term' 
in a set of equations of motion of the De Witt/ Brehme p] form. Direct 
mathematical comparison with the equations of motion in pQ is highly non- 
trivial since the equations obtained in pQ have involved the removal of an 
infinite term, while (I3.74p does not involve such a procedure. Hence we say 
that ( 13.741) fits the pattern of the equations obtained in pQ. 

4 Discussion 

Perhaps the most significant aspect of the work presented in this paper is the 
fact that the equations of motion ( 13.741) of a small charged black hole moving 
in an external electromagnetic and gravitational field have been derived from 
the vacuum Einstein-Maxwell field equations without encountering infinities. 
It is also important to note that the external fields have been introduced as 
a solution of the Einstein-Maxwell vacuum field equations. 

The second term on the right hand side of ( 13.741) is the electromagnetic 
radiation reaction 4-force. The third term on the right hand side of (13. 741) 
can be written in terms of the electromagnetic energy tensor Eij{u) of the 
background space-time, calculated on r = in coordinates {X 1 ). Using the 
background field equations on r = 0, i2y = 2Eij, this term can be written 
in terms of the background Ricci tensor components as (2/3)e 2 /if R k j . 
This is the 4-force due to the external field and it differs from that derived 
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by Hobbs [35] by a factor of 2. The external electromagnetic field in [33] 
is not required to satisfy the Einstein-Maxwell vacuum field equations and 
so the background Ricci tensor appearing in [35] does not have the external 
electromagnetic field as source. Thus in effect half of the contribution to the 
external 4-force is neglected. 

The 'tail term' in the form (I3.75P vanishes if the external electromagnetic 
field vanishes. It involves three arbitrary functions (of u) K, U, V with the 
latter two obtained from U, V, X, Y as indicated following (13. 73ft . The math- 
ematical origin of these functions as functions of integration arising in the 
determination of K-i, h and m-i, required for the perturbed 4-potential, and 
of cio and bo required for the perturbed metric tensor, is clear from this work. 
In addition the tail term indirectly involves the two vector fields c* and d l 
orthogonal to v % . It is a topic for further study to understand physically 
the role of these arbitrary functions and of these unique space-like vectors. 
Had the tail term not been so explicit we could not identify these arbitrary 
functions and space-like vectors for further consideration. By comparison 
the tail term obtained in pQ is an integral whose integrand is expressed in 
terms of functions appearing in the Hadamard form of the Green function of 
the vector wave equation. Very little is known explicitly about the form of 
such functions. The De Witt-Brehme tail term has proved quite intractable 
to analyse due to the fact that the Green function is not known in closed 
form. This makes a comparison with our result extremely difficult. Such a 
study would surely merit a paper independently of ours. 

The original calculations of Dirac [2] were carried out in Minkowskian 
space-time. Our results do not specialise to those of Dirac since our back- 
ground space—time is an Einstein-Maxwell space-time and if it is flat then the 
external electromagnetic field vanishes. In this case f!3.74p gives maj = 2 
and when this is substituted into the surviving radiation reaction 4-force the 
latter is absorbed into the 03-error. Hence we obtain geodesic motion (and 
in particular no "run-away" motion) at this level of approximation in this 
theory. 
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A Useful Formulas from Minkowskian Geom- 
etry 

We record here some useful equations which hold exactly in Minkowskian 
space—time and are also important in the neighbourhood of the world line 
r = in the background space-time in this paper. The transformation (I2.24p 
can in principle be inverted giving (x, y, r, u) as functions of (X 4 ) which means 
that we can consider (x, y, r, u) as scalar fields on Minkowskian space-time. 
Their derivatives with respect to (X % ), denoted by a comma, are obtained 
by first differentiating (12.241) to arrive at 



8] = (V - r h k l ) u tj + k l r tj + r ~g^ x J +r ~g^ Vj ■ ( A - l ) 



Multiplying this by hi gives immediately 

kj = Uj . (A.2) 

Now multiplying ( 1A.1I) by Vj yields 

r ,j = v 3 - (! -rh )kj . (A.3) 

Differentiating ( 12. 24ft with respect to X J and using ( 1A.2I) and ( 1A.3I) results 
in the alternative form for ( 1A.1I) : 

k\j = - (6) - v l kj - k l Vj + (1 - r h ) k % kj) , (A.4) 

with ho = dik 1 . From this it follows that 

-2L = jfe* .y = , — = k* j v j = -h k* . (A.5) 
or du 

Writing ho = <9(log Po)/du for some function Pq independent of r we see 
from the latter two equations that k l = P _1 ( l with ( l null and independent 
of r and u. Thus Q can be parametrized by two real parameters. We have 
chosen the parameters x, y as given in (I2.23f) and then Pq = Q v l in ( 12.201) is a 
consequence of this choice. A more complete discussion of this construction 
can be found in [33] and [34] . 

We can now simplify flA.10 to read 

dk l dk 1 

5 l - = v % k j + k l v.-, - k l kj + r — — x.j + r ^— y j . (A.6) 

ox ay 
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Multiplying (IA.6I) by dki/dx and also by dki/dy provides the remaining 
equations needed: 

x ,j = ~l P o ^ > Vj = ~Po ^ • (A.7) 

Substituting f]A.7|) into f1A.6[) results in the Minkowskian metric tensor being 
written in the form 

^ = " P o f ^ ^ + ^ f^l + fcJ vj + k3 «' - ki kJ > (A-8) 
\ ax ax a?/ ay / 

and much use is made of this relation in the calculations behind this paper. 
Equivalently the Minkowskian line-element takes the form 

ds 2 = r]ij dX i dX j = -r 2 P^ 2 {dx 2 + dy 2 ) + 2 du dr + (1 - 2 h r) du 2 . (A.9) 

This suggests we introduce basis 1-forms (defining a half null basis) 

tf 1 = r P 1 dx , $ 2 = rP 1 dy , $ 3 = dr + ~(1 -2 h r) du , tf 4 = du . (A.10) 

Using ( 1A.2I) . ( 1A.3I) and ( 1A.7I) we can express these in terms of the rectangular 
Cartesian coordinates and time as 

Bk- 

V 1 = -P o -^dX i = -0 1 , (A.ll) 

dh- 

§2 = - p ^ dXl = - § ^ (A- 12) 

tf 3 = (y i -^k i )dX i = &i, (A. 13) 

tf 4 = kidX l = d z . (A. 14) 

It is helpful to have available the second partial derivatives of k % with re- 
spect to x and y expressed on the basis consisting of the vectors v\ k\ dk l jdx 
and dk % jdy. These formulas are: 

, (A.15) 
, (A.16) 
(A.17) 



d 2 k l 

dx 2 


= Po 2 (v l - 


-V) 


d n 


Po) 


dt 
dx 


d n 

+ ^ 


;Po) 


dk l 
dy 


d 2 k { 
dy 2 


= Po 2 ^- 


■ k l ) 




Po) 


dk l 
dx 




;Po) 


dk l 
dy 


d 2 k l 
dxdy 


8 n 
^ -dy- {l ° g 


Po) 


dk % d 
dx dx 


(log 


Po) 


d& 
dy 
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B Perturbed Field Equations 



Writing, for the perturbed space-time described in section 3, M. a = F ab \ b and 
£ ab = Rab — 2 E ab the leading terms and the errors we tolerate in these expres- 
sions (and thus the extent to which we satisfy the vacuum Einstein-Maxwell 
field equations) are given here. If the equations of motion are required with 
greater accuracy then the field equations have to be solved with greater ac- 
curacy resulting in smaller errors in the coefficients of these powers of r. 

M A = 3 x r- 2 + 2 x r~ l + 1 + 0(r) , (A = 1,2), (B.l) 

M 3 = 3 x r- 1 + Ox + 0(r) , (B.2) 

M A = 2 x r + 0(r 2 ) , (B.3) 

and 



£ 



£aa 


= 3 x r~ 4 + 3 x r~ 


" 2 + 2 x r- 1 + O(r ) , 


(B.4) 


11 — £22 


= 3 x r~ 2 + 2 x r~ 


~ 1 + 1 + 0(r) , 


(B.5) 


£12 


= 3 x r -2 + 2 x r~ 


- 1 + 1 + 0(r) , 


(B.6) 


£A3 


= 3 x r -2 + Oi x 


- 1 + 1 + 0(r) , 


(B.7) 


£aa 


= 2 x r~ 3 + 2 x r~ 


- 2 + 2 x r- 1 + O(r ) , 


(B.8) 


£33 


= Oi + 0(r) , 




(B.9) 


£34 


= 4 x r~ 4 + 3 x r" 


" 3 + 2 x r~ 2 






+0 2 x r- 1 + 0(r°) 




(B.10) 


£44 


= 3 x r~ 3 + 3 x r" 


" 2 + 0(r- 1 ) . 


(B.H) 



The coefficients of the various powers of r in these expressions are calculated 
by substituting into A4 a and £ ab the metric tensor, given via the line-element 
(12.11) with the expansions ( 13.ip - 113.6p . and the potential 1-form (12.141) . with 
the expansions (I377j) - (13T3]) . Writing (!B~ll - (lB~Tll in the form 

M a = J2(n)M a r n , £ ab = J2(nA b r n , (B.12) 

the coefficients ( n )Ai a and ( n )£ a & required to establish ( IB. If) ( TB. 1 If) are: 



9 Jp 2 fdM n dL 



. )r +0 3 , (B.13) 
(-i)-M 1 = 3 , (B.14) 



(o) 



M 1 = QmL 3 + 2 = 1 , (B.15) 
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(-2) 



M 2 



(o)-M 2 

(o)-M 3 
(i) 



M 4 



Writing 



we find that 



{-a)£aa 
(-3)£aa 



AA 



(-2)^11 — (-2)^22 



(-1)^11 — (-1)^22 



... p 2 fdM dL 



d_ 

_dK-i 
dy 



\ dx dy 
+ 3 , 



2e 2 M 2 + eb^pQ 



03, 

-6mM 3 + 2 = Oi , 
2 , 

-4mPo 2 F ij k i v j + 02 = 0! , 

-2P - 2 £Ti + AP 2 (a_iL 2 + S_iM 2 - eq 2 ) 

(dL 2 dM 2 \ „ 
+2 ( -z 1 + I +O2 ■ 



<9x 



- / <9 2 <9 2 
A = P 2 (^ + 



(9a; 2 <9y 2 



3 , 
3 , 



(B.16) 

(B.17) 
(B.18) 
(B.19) 
(B.20) 



(B.21) 



(B.22) 



(B.23) 
(B.24) 



{ - 2) £aa = 2AlogP -4P 4 



'dM dL \ (dM 2 dL 2 



dx 



dy 



dx 



dy 



_ 3 P* (^a.,) + | ( p-6_,)) 

+4eKi + 4Ki X_i -8e(a_iL 2 + 6-iM 2 ) 
-2 c + 12 e 2 g 2 - ^P 2 (a 2 _ 1 + b 2 _ J + 3 , 

-4ci-8 " g ° +8eif 2 -32mg 2 
cm 

- 4 ^(cl ( ^ 2ao) + c| ( ^ o) ) +02 ' 
^-^-4e« 2 --P (a_ 1 + M 



-S, 2 /^/. 2 . Mf) 8P 2 L 2 5A 



+0 3 , 



dx 
dK 



+ 8P 2 M 2 



16mP 2 (Lj-M 2 )-8P 2 L 2 ^ + 8P 2 M 2 ^ 



(B.25) 



(B.26) 



dy 

(B.27) 
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+16m« 2 -2^ + 2^ + 2 , (B.28) 
ox ay 

-8c P 2 (L 2 - Ml) - 8 PfL 2 ^ + 8 P 2 M 2 ^ 

-12 a 2 c -3-^ + 3-^ + O 1 =O 1 , (B.29) 
1 (9a_i 1 <%i 9 A 1 a 

-8 e 2 P 2 L 2 M 2 - 4 P 2 L 2 - 4 P 2 M 2 

+0 3 , (B.30) 

16 m P*L 2 M 2 -A PlL 2 - 4 P 2 M 2 

+8m/? 2 -— — h 2 , (B.31) 

ay ax 

-8c„ P % M, - 4 P|M 2 ^1 - 4 P 2 L 2 ^1 

- 6 A--^-|| + 0l = 0,, (B.32) 

4M 2 P =(f -^) + 4ei 2 P + Va- 

+4K_ 1 L 2 Po + 3 , (B.33) 

6eP L 3 + O 2 = d , (B.34) 

-4 ^ L 2 P + 2 P 3 ( J^o" 2 «2) + |^o~ 2 &)) 

= 1 , (B.35) 

9 A;* 

4meP ^^ + O 3 = O2 , (B.36) 
2 e c L 2 P + 2 e P)^ 1 + ^ P^a-i c 

-5<(^{|^ ,4 -''-s^ L ')}) 

-P _1 a-i A log P + 2 = 2 , (B.37) 
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(-i)£i4 = 2 , (B.38) 

= -^ 2 n 3 (^-^)+4eM 2 P + P - 1 ^ 

+4X_ 1 M 2J P + 3 , (B.39) 

( _i)£23 = 6ePoM 3 + 2 = Oi , (B.40) 

( o)^ 2 3 = -4K 1 M 2 P + 2P 3 ^-^(P - 2 d 2 ) + ^(P - 2 ^)^ 

= Oi , (B.41) 

( _3)£ 24 = 4 m e? F {J -^+O 3 = O 2 , (B.42) 

(-2)^24 = 2ec M 2 Po + 2ePo^ + -P ' 1 L 1 co 

-P _1 6-i A log P + 2 = 2 , (B.43) 

(-1)^24 = 2 , (B.44) 

( o)^33 = -8P 2 (Z| + M 2 2 ) + 12g 2 + O 3 , (B.45) 

(-4)^34 = 4 , (B.46) 

(-3)^34 = 3 , (B.47) 

{ _ 2) S U = 2e^ 1 + i^^(P - 2 a_0 + | ; (Po- 2 S-i)) 

+0 2 , (B.48) 
<9P 

(-1)^34 = Ci+2P - 1 ^ + 4ei 2 + 8mg 2 

+ ^ ( J^ "^ + ^ ( ^ o) ) + ° 2 ' (B ' 49) 

(-4)^44 = O s , (B.50) 

M) ^44 = 2mP 2 (i(P 2 L 1 ) + |(P l 1 )) 

an 

-4e 2 Po-^-2A/_ 1 + 03 (B.51) 
aw 

( _ 2) ^44 = T\ + T 2 + T 3 + 3 . (B.52) 
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In the final equation here the terms T 1; T 2 , T 3 are given in the text by (13.441) . 
(I3.55P and (13.571) respectively. 
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